Collective excitations of graphene excitons being in the Bose-Einstein condensate state 
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Bose-Einstein condensation of the excitons in graphene is considered. We suggested the model spinor 
equation for neutral particles with short range interaction described the microscopic graphene excitons 
dynamic. Using this equation we derived quantum hydrodynamic equations for description of collective 
properties of excitons in graphene, particularly for the case when excitons being in the Bose-Einstein 
condensate (BEC) state. The dispersion of collective excitations in graphene excitons BEC is studied. 
We shown that frequency of collective excitations is proportional to the square root of the wave vector 
module: cu ~ yk. 



Bose-Einstein condensation (BEC) of ultracold atomic 
gases 0]] and excitons in semiconductors [0]- J@] and 
monolayer of carbon atoms (graphene) 01, J8|]7 10] 
are very interesting and promising physical systems for 
present day. In this paper we present model for excitons 
in graphene and study they properties in the case one's is 
in the BEC state. In connection with it we make some no- 
tice about kinematic and dynamic properties of particles in 
various fundamental models. 

Quantum mechanics and special relatively gives us a new 
way of the kinematic properties description in comparison 
with the classic physics. However, the methods of inter- 
action account remains the same as in classic physics. For 
example, in the Schrodinger equation the Coulomb interac- 
tion is considered as interaction of point like particles IToll . 
Some times the Coulomb interaction is also used in the 
many-particle generalization of the Dirac's equation fiUl . 

In graphene we have deal with the electrons and holes 
conductivity. They properties crucially depend on fields 
caused by the lattice. In Ref. Il2ll were found that mo- 
tion of carriers in graphene might be described by Dirac's 
like equation. Therefore, we met a new type of kinematic 
properties when we study graphene carriers. In Ref. IU3I1 
the Coulomb interaction were added in Hamiltonian for 
graphene electrons. In this paper we consider the short 
range interaction (SRI) between excitons in graphene and 
we suppose that the graphene excitons has the same kine- 
matic properties as the graphene carriers. 

Dirac's equation describe the relativistic motion of elec- 
tron in external electromagnetic field. This equation so- 
lution for energy eigenvalues give us solutions with posi- 
tive and negative energy separated by energy gap. Solu- 
tion with the negative energies corresponds to the antipar- 
ticle states; i.e. to the positron states. When we have deal 
with the semiconductors we work with the electron's and 
hole's. Holes are the quasi-particles corresponds to the mo- 
tion "ionized" states of atoms in semiconductor. To move- 
ment of hole corresponds the movement of bound electrons 
in opposite direction. Electrons of conductivity moving in 
the space between atoms as quasi-free particles. To the hole 
motion correspond exchange of electron from one bound 
state to another between neighbor atoms. Thus, mecha- 



nisms of described two types of motion are different. This 
leads to the different effective mass of electrons and holes, 
and it's difference from mass of the free electron. 

In the literature there is the analogy between behavior 
of relativistic electrons (the picture of electron-positron 
states) and electrons in semiconductors (the picture of con- 
ductivity electrons and holes). Electrons and positrons has 
the same mechanism of motion. Difference in the mech- 
anisms of motion and in the effective masses of the elec- 
tron's and hole's lead to the fact that using of Dirac's equa- 
tion for description of the electrons with the positive and 
negative (hole) energy level is not always suitable for semi- 
conductors. 

Nevertheless, massless Dirac's like spinor equation has 
been used for the description of the conductivity electrons 
in graphene 



(1) 

where ip is the A r -particle wave function which depends on 
2N coordinate, because graphene is the two-dimensional 
(2D) structure, v F is the Fermi velocity of conductivity 
electrons in graphene, electric charge, a a is the spin 
1/2 Pauli matrixes, in this equation index a attain to x 



and y, D { 
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the scalar and vector potential of external electromagnetic 
field, Gij = 1 /r jj is the Green function of Coulomb inter- 
action. The same equation might be used for description of 
the holes in graphene. 

Electrons and holes are the charge carriers in graphene. 
They can form bound states - excitons. Typical velocity of 
electrons and holes in graphene is the Fermi velocity. The 
mechanisms of motion of electrons and holes in semicon- 
ductors are different, but the typical velocity of they motion 
in graphene is the same. Thus, we can suppose that center 
of mass of excitons in graphene move also with the Fermi 
velocity. Important difference of excitons from electrons 
and holes it is the fact that spin of exciton equal to inte- 
ger number of Plank constant: or h. Excition are Bose 
particles, whereas electrons and holes are Fermi particles. 

As the atoms in quantum gases as the excitons in semi- 
conductors attach a lot attention in the connection with the 



2 



Bose-Einstein condensation. For the first time BEC was 
realized in vapors alkali atoms in 1995. In atomic gases 
the BEC state reached at hundreds of nanokelvin, whereas 
in exciton systems BEC state might be realized at several 
kelvins Jl, Q, Q, d. 

We suppose it is possible to realize the BEC of excitons 
in graphene and suggest model Dirac's like equation for the 
neutral particles with short range interaction for description 
of excitons in graphene. This equation is an analog of the 
equation suggested in Ref.s II 1211 . II 1 311 . Here we present 
some notice about using of Dirac's like equation for Bose 
particles description. Used in Ref. Ill3ll equation is linear 
on momentum operator as Dirac equation. Instead of scalar 
product of Dirac's matrixes on four momentum this equa- 
tion contain the scalar product of Pauli matrixes [<j x , a y ] on 
momentum operator p = p j. Thus, the model Dirac's 
like equation used in Ref. 111311 might be considered as an 
analog of spinor Pauli equation containing another Dirac's 
like kinematic. Consequently, such model equation as the 
Pauli equation might be used for the description as Bose as 
Fermi particles. 

In this letter we found the general equation for the col- 
lective motion of the system of neutral spin-1 particles with 
SRI whose kinematic properties is analogous to graphene 
carriers. From this general equation we derive equation for 
BEC in described particles system. For derivation of equa- 
tions described collective motion from the many-particle 
Dirac's like equation we use the method of many-particle 
quantum hydrodynamics (MPQHD). This method was sug- 
gested in 1999 by L. S. Kuz'menkov and S. G. Maksimov 
1 17y. Further development of the MPQHD was made in 
Ref.s lfl8ll- i24fl fo r various system of particles: quantum 
plasma I17ll. Il8ll. rf2~lll . particularly for spining parti- 
cles 11181 . 1211 1. 1221 1: BEC and ultracold fermions of neu- 



tral atoms 119D; charged and neutral particles with elec 
trie dipole moment [20], particularly electrically polarized 
BEC [23]; and electrons in graphene [24]. M. Marklund 
and G. Brodin suggested another way of derivation of QHD 
equations for spining particles lE5ll . JU. 

The QHD equations as a classic hydrodynamics are very 
useful for collective excitation studying Il20ll . I2lll . 112311 - 
1251. J28]. In this paper we consider the spectrum of col- 
lective excitations in BEC of excitons in graphene. 

In this letter we do not present the detail of derivation 
of QHD equation for excitons BEC in graphene (GEBEC). 
The general scheme of QHD equations derivation is pre- 
sented in Ref.s d, H3], and I24ll . 

We consider BEC in 2D system of particles. BEC can- 
not take place in a purely 2D system at finite temperatures. 
Consequently this term may be used for polaritons under 
very special circumstances (see 1E9I1 and ref.s where). De- 
tailed discussion of the superfluidity phenomena and BEC 
in 2D system is presented in Ref. [30]. Influence of dipo- 
lar nature of excitons on one collective properties is also 
studied where. 

For exciton-polariton BEC description the Gross- 



Pitaevskii equation is used 113 ill as for exciton-polaritons 
dynamic in non-condensed states lf32h . Gross-Pitaevskii 
equation is usually used for studying of BEC in atomic 
vapors [1]. Gross-Pitaevskii equation is the non-linear 
Schrodinger equation. This equation has form of the 
Schrodinger equation for one particle in external field, 
but also contains the non-linear term proportional to the 
third degree of the wave function. Gross-Pitaevskii equa- 
tion might be present in the form of two hydrodynamic 
equations: continuity equation and Euler equation yj]. 
Method of QHD allows to derive the Gross-Pitaevskii 
equation and it's nonlocal generalization from many- 
particle Schrodinger equation Ill9ll . 

We present a basic equation here for evolution descrip- 
tion of excitons in graphene. Using this equation we derive 
and present below the equations for description of exci- 
tons collective motion and especially for the dynamics of 
GEBEC. 
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this equation differ from (fl]l by the form of spin matrix and 
the form of interaction. Equation (O contains following 
quantities: wave function ip = ip(R,t), R is the whole 
particles coordinates R = [ri,...,rj, ...,rjv], r$ = [x,,^], 
sf are the spin-1 matrixes for i-th particle, pf = —iftV 
is the momentum operator, Vi iext is the potential of exter- 
nal field, Uij is the SRI potential describing the interaction 
between excitons in graphene. We notice that at the same 
time some particle might interact with several particles, by 
means SRI potential . For this statement illustration we 
refer to the liquid where molecules is neutral and interacts 
with the several neighbor molecules. We consider spin-1 
particles and spin operators are 3 x 3 matrixes 
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the commutation relation for spin-1 matrixes is 

[s«,s^=ikie aM sl. (3) 

Equation ([2} describes the 2D motion of excitons, thus 
Hamiltonian contains operators Sf, Sf, p x t and p\ only, but 
at using commutation relation ([3]l during equations deriva- 
tion the -operator is also appearing. 

Graphene is the 2D structure and electrons of graphene 
are located in the plane. As we describe above in 2D case 
the electrons has two coordinate x and y, but spin of elec- 
trons can be directed in all direction, particularly, in z axes 
direction, perpendicular to the graphene plane. This fact is 
accounted by formula (O. 
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The first equation of the QHD equations set is the conti- 
nuity equation 

d t n(r,t) + v F VS(r,t) = 0. (4) 

The definition of probability density of conduction electron 
system in physical space is 

n(r,t) = J2 [ dRY,8(T- ri W(R,t)i/i(R,t) (5) 

where dR = n»Li dr p . The continuity equation appears 
at differentiation of the concentration Q with respect to 
time and using equation ® corresponding to the QHD 
method. 

The quantity n(r, t) can be considered as 2D excitons 
concentration. The time evolution of concentration The 
quantity S(r, t) describe the spin density of the system of 
particles. Coordinate vector r has only two component. 
Consequently, equation © contains two component of spin 
density vector S, these are S x and S y . Our next step in 
construction of the model of collective motion is obtaining 
of equation of spin evolution. For this aim we differenti- 
ate quantity S(r, t) with respect to time and use equation 
©. Because we known the third component of spin den- 
sity vector we can study the evolution of whole component 
of this vector. Therefore, we derive evolution equation for 
S = [S x , S y , S z ]. On this way we have equation of spin 
evolution: 



d t S a (r,t)+v F 
1 



n(r,t) 



(6) 



Here new physical quantity is appeared: JJjf (r, t). This is 
the tensor of spin current. 



n 2 (r,t) 



+f £ ^(r,t) 

a \ n z {r,t 



J»(r,t)JP(r,t) 2 d^d ^/n(r^) 



dr>(d?U(r, r'))j 2 a (r, r', t)-v F S a (r, t)d V ext (r, t). 

(7) 

Equation (O contains interaction and a new quantity 
J a (r, t), J a (r, t) is the probability current. 

d t J a (r,t) +v F d J^ a = -v F f dr'(^[/(r,r'))x 

x nj(r, r', t) - v F n(r, t)d a V ext (r, t). (8) 

The equations (0 and (HJ contains two-particle function. 

These two-particle function appear in the terms contains 
SRI. In Ref. lfl9ll was developed the method of calcula- 
tion of the this kind quantities. Using two condition: par- 
ticles interact via SRI and particles being in the BEC state, 



we obtain closed QHD equations set. We do not stop here 
on the technic of many-particles function calculation. For 
neutral particles being in the BEC state this technic was 
described in the Ref. 11911 . Here we present the result- 
ing equation found at the condition that excitons are in the 
BEC state. 

d t J a + VF&jfr = v F T 2D nd a n - v F nd a V ext (9) 
and 



,p , v F /5 Q 5^Jn 1 



dt Jm P + fd 



-v F T 2D df3(S a n) 



+ 



h 



n- 



n 



v F 



S a d V r ,, 



(10) 

In these formulas we do not present the arguments of func- 
tions. 



T 



7T 



dr(r) 



2 dU(r) 



dr 



(11) 



this quantity is the analog of the interaction constant in 
Gross-Pitaevskii equation 0]], fl9Tl . In the BEC theory the 
SRI interpret via scattering, and T = —4irh 2 a/m, where 
a is the scattering length in first Born approximation yj]. 
In Ref. 12311 were presented that at derivation of QHD 
equations we do not consider the scattering approximation 
of interaction and connection of T and scattering length a 
presented here for handly comparison. 

We can analyze the linear dynamics of elemental excita- 
tions in the polarized BEC using the QHD equations 0, 
©, (O and ( fTOl ). The system is placed in an external mag- 
netic field B = B e z . The values of concentration n and 
polarization S || B for the system in an equilibrium state 
are constant and uniform and its velocity field J Q (r, t) and 
tensor J Q ^(r, t) values are zero. 

We consider the small perturbation of equilibrium state 

n = n Q + 5n, S a = S« + 5S a , 

J a = + 5J a , J af3 = + 5J al3 . (12) 

Substituting these relations into system of equations 0, 
©, © and ( fTUl) and neglecting nonlinear terms, we ob- 
tain a system of linear homogeneous equations in partial 
derivatives with constant coefficients. Passing to the fol- 
lowing representation for small perturbations 5f 

Sf = f(uj, k)exp(—iujt + ikr) 

yields the homogeneous system of algebraic equations. 
The spin density strength is assumed to have a nonzero 
value. Expressing all the quantities entering the system 
of equations in terms of the spin density, we come to the 
equation 

A a0 (u,k)SP{u>,k) = 0, (13) 

Excitations exist in the case when matrix equation dT3l) has 
nontrivial solutions. The condition of nontrivial solution 
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existence for homogeneous linear algebraical equation set 
is the determinant of this equation set must be equal to zero: 

det || A Q/3 (cj,k) ||= 0. (14) 

Solving this equation with respect to u 2 we obtain a fol- 
lowing results. 

Equation (fl4l give us only one wave solution, dispersion 
of this wave is obtained in the form: 



v F n 



T | k 



hV2 



(15) 



where k 2 = k 2 + k 2 . 

Solution (Q3J does not depend on the sign of SRI, and 
contain module of the SRI constant T. 

For comparison we present here the dispersion of the Bo- 
goliubov's mode which is the eigenwave in usual BEC: 



2 ^ ; 4 , n , 2 

u = ——k +g—k , 

4m- 2 m 



(16) 



in the long wavelength limit from formula ( fl"6l) we have 
uj 2 = gn k 2 /m. The quantity g is the interaction con- 
stant in Gross-Pitaevskii equation [1] and g = — T. This 
solution exist for repulsive interaction (g > 0, a > 0) only. 

Thus, instead of sound like spectrum uj ~ k for BEC in 
atoms vapors for the GEBEC we have lo ~ Vk. 

In connection with the studying developments of the 
BEC in excitons system in semiconducturs and wide us- 
ing of graphene in semiconductor geterostructures we sug- 
gest the model for description of GEBEC collective prop- 
erties. One of the fundamental properties of many-particle 
systems is the spectrum of collective excitation. Thus, we 
considered the collective excitation spectrum for GEBEC. 

Using the idea that excitons in graphene has an analo- 
gous kinematic properties as the graphene carriers we for- 
mulated many-particle Dirac's like equation (or linear on 
momentum pauli like equation) for spin-1 neutral particles 
with SRI, describing graphene excitons dynamics. Starting 
from this equation, using method of many-particle QHD 
we derived the system of QHD equations: particles number 
balance equation (continuity equation), spin balance equa- 
tion, current evolution equation (Euler equation) and equa- 
tion evolution of spin current. We considered graphene ex- 
citons being in the BEC state and studied the dispersion of 
collective excitations where. We found that derived model 
leads to existence of one type of collective excitations in 
GEBEC. Dispersion dependence of these excitations dif- 
fer from Bogoliubov's spectrum, which suitable for BEC 
of atomic vapors or liquid helium. In the case GEBEC 
the frequency of collective excitations is proportional to 
\fk. Thus, obtained dependence differ from sound wave 
and long wave length limit of Bogoliubov's spectrum. 

The author thanks Professor L. S. Kuz'menkov for fruit- 
ful discussions. 
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